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A model for studying atomtronic devices and circuits based on finite temperature Bose-condensed 
gases is presented. The approach involves numerically solving equations of motion for atomic pop¬ 
ulations and coherences, derived using the Bose-Hubbard Hamiltonian and the Heisenberg picture. 
The resulting cluster expansion is truncated at a level giving balance between physics rigor and 
numerical demand mitigation. This approach allows parametric studies involving time scales that 
cover both the rapid population dynamics relevant to non-equilibrium state evolution, as well as 
the much longer time durations typical for reaching steady-state device operation. The model is 
demonstrated by studying the evolution of a Bose-condensed gas in the presence of atom injection 
and extraction in a double-well potential. In this configuration phase-locking between condensates 
in each well of the potential is readily observed, and its influence on the evolution of the system is 
studied. 


I. INTRODUCTION 

In many respects atomtronic devices and circuits mm 
can be understood using the same framework applica¬ 
ble to their electronic counterparts. Analytical tools 
such as Kirchhoff’s voltage and current laws, for exam¬ 
ple, are simply circuit-relevant statements about energy 
and particle conservation, which apply equally well to 
electronic and atomtronic systems. Operating in the ul¬ 
tracold regime, two aspects of atomtronic circuits make 
them both interesting and very different from classical 
electronics. First, quantum coherence can be long lived 
compared with other circuit timescales of interest. Sec¬ 
ond, atomtronic circuits are typically decoupled from a 
thermal reservoir [3]. An active circuit is by definition a 
non-thermal-equilibrium dynamical system and any use¬ 
ful circuit will generate heat for both fundamental and 
technical reasons. In general, this means that the tem¬ 
perature will vary from place-to-place across the circuit. 
Thermal coupling to the environment enables an elec¬ 
tronic circuit to reach thermal steady-state (not thermal 
equilibrium) and finite device temperature. However, 
heat has no thermal bath to flow to in an atomtronic 
circuit. While doubling its temperature often leads to 
permanent failure of a nominal room-temperature elec¬ 
tronic device, the same is not true of an atomtronic de¬ 
vice. Yet such significant changes in temperature can 
certainly have an impact on the quantum character of 
devices operating in the ultracold regime. The interplay 
between the thermal-statistical and quantum aspects of 
atomtronic circuits at the physical level is profoundly re¬ 
lated to circuit functionality through entropy - physical 
entropy on the one hand and information entropy on the 
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other. For example, if a circuit function corresponds to 
the lowering of information entropy, as is often the case, 
this must correspond to a lowering of physical entropy in 
some region of the circuit, which demands that there is 
a corresponding increase in entropy in another region of 
the circuit if the second law of thermodynamics is to be 
upheld. 

Electronic circuits are challenging from an analytical 
standpoint and atomtronic circuits are dramatically more 
so. Only the most elementary electronic circuits can 
be treated from first principles; complex circuits are de¬ 
signed using heuristic principles and analyzed using pa¬ 
rameterized device models and highly optimized numer¬ 
ical solvers. Our purpose in this work is to make some 
headway along these same lines in the atomtronics do¬ 
main. We motivate our objective with a few simple ques¬ 
tions regarding a conceptually simple system. As it turns 
out, the questions are very difficult to answer. 

Consider the double-potential-well system of Fig. [lja). 
The right-hand well is initially devoid of atoms while the 
left-hand well holds a Bose-Einstein condensate (BEC) 
at, or very near to, absolute zero. The BEC is character¬ 
ized by a chemical potential, while the empty right- 
hand well has a chemical potential equal to its ground 
state energy. For this discussion we take fx l to be large 
compared to the single-particle energy level spacing near 
the ground state of the well, but considerably less than 
the barrier height between the wells. Our first question 
is: “Does this (isolated) system come to thermal equi¬ 
librium, and, if so, by what route does it do so?” By 
thermal equilibrium we mean the system is described by 
a single, time-independent temperature, T, and chemical 
potential, /r, as illustrated in Fig. 0b). We note that the 
equilibrium state must have finite temperature to main¬ 
tain fixed total internal energy, even if the initial left- 
hand BEC is at zero temperature. Thus, by whatever 
route the system arrives at equilibrium, heating must 



2 



C L3 >» 

U) 

E R2 ® 

c tU 



1 ”L3 >> 

O) 
E R2 o 

c UJ 


FIG. 1. Schematic representation of a BEC heating as it 
comes to thermal equilibrium in a double-well potential, (a) 
Initially, the state of the system is highly non-equilibrium as 
demonstrated by the empty states in the population distri¬ 
bution at energies corresponding to levels in the right-hand 
well, (b) In thermal equilibrium the population has been re¬ 
distributed via elastic collisions between atoms because pop¬ 
ulation of lower lying states in the right-hand well are always 
accompanied by population of higher lying states. The net 
result is heating of the initial distribution as depicted by the 
change in color from (a) (blue, colder) to (b) (red, hotter). 


take place. The energy associated with the temperature 
rise comes at the expense of chemical potential, which 
necessarily decreases from its initial value in the left- 
hand well. Supposing that parameters are chosen such 
that thermal equilibrium supports a condensate in both 
wells, then their chemical potentials must be equal in 
thermal equilibrium. This implies that their respective 
wavefunctions have a well-defined relative phase which, 
ignoring fluctuations, is constant in time. Our second 
question is: “Can the two condensates phase-lock, and if 
so, under what conditions does this occur?” By “phase- 
lock” we mean that if the relative phase is perturbed it 


subsequently returns to a fixed equilibrium value. Phase¬ 
locking is a ubiquitous phenomenon in nature (see [4] for 
a number of examples), and more to our point, it is an im¬ 
portant functional primitive among electronic circuits - it 
lies at the heart of every modern radio receiver, for exam¬ 
ple. In the circuit context, questions regarding thermal 
equilibrium can similarly be asked regarding dynamical 
equilibrium. 

To make an attempt at answering both questions in 
the circuit context, we developed a model of an atom- 
tronic circuit operating with finite temperature BECs 
in the presence of atomic current flows, which are de¬ 
termined by atom injection and extraction processes, 
as well as collisions. The model draws upon a wealth 
of theoretical developments involving weakly interacting 
Bose gases [SHE]. Past studies of phenomena occurring 
in multi-well potentials include atom interferometry [8], 
the bosonic Josephson effect [5], dynamical phase-locking 
of BECs [10], and transistor-like behavior of ultracold 
atoms m- They typically involve solutions of the time- 
independent Schrodinger equation, or the more general 
von Neumann equation for the Bose-Hubbard Hamilto¬ 
nian [12] . and treat collision effects at the level of the 
Boltzmann equation [3 Eg. In these studies, detailed 
analyses were performed, sometimes under restrictive ex¬ 
perimental conditions, such as for a microcanonical en¬ 
semble with a fixed total number of atoms to isolate the 
phenomenon of interest in order to understand its under¬ 
lying physics. The methods are effective at extracting in¬ 
trinsic and extrinsic parameters that characterize the gas, 
such as damping rates for various collective excitations. 
These parameters provide good anchors for our model, 
which is more phenomenologically based for two reasons. 
First, the number of atoms in an atomtronic circuit is 
not known a priori and is likely not conserved. There¬ 
fore, one has to solve the more general time-dependent 
problem accounting for the influence of atom injection 
and extraction on BEC dynamics. Second, the neces¬ 
sary computations cannot be so demanding as to prevent 
timely and comprehensive parametric studies. 

Our starting point remains the Bose-Hubbard Hamil¬ 
tonian. However, we choose to work in the Heisenberg 
picture and use a cluster expansion approach to circum¬ 
vent the complication of a rapidly increasing configura¬ 
tion space with an increasing number of atoms. The 
approach involves tracking the time evolution of atomic 
populations in the stationary states of the confining po¬ 
tential. Many-body equations of motion for atomic popu¬ 
lations and coherences are derived. They contain hopping 
effects, energy renormalizations, and coherences within 
the context of a mean-field approximation. Collisions 
giving rise to relaxation and dephasing are modeled phe¬ 
nomenologically using an effective relaxation rate approx¬ 
imation, with the effective relaxation rates obtained from 
experiment or quantum kinetic calculations |7j, for exam¬ 
ple. This type of treatment is common in quantum elec¬ 
tronics M- Beyond the investigation of phase-locking, 
our interest is to understand the extent to which the con- 
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FIG. 2. Confining potential (solid, black) for the double-well 
geometry being studied. The parabolas (dashed, black) rep¬ 
resent the potentials for the uncoupled wells. Also plotted are 
the eigenfunctions (solid, blue and red) of each well displaced 
vertically according to energy. 


siderable body of work in many-body semiconductor laser 
theory mm can provide tools for designing and ana¬ 
lyzing atomtronic devices and circuits. 

Within our model a trade-off is made between rigor and 
mitigation of numerical demands. The result is an ap¬ 
proach allowing parametric studies involving timescales 
spanning over three orders of magnitude. With this 
dynamical range rapid population dynamics relevant to 
non-equilibrium effects can be modeled, as well as the 
long time durations typical of steady-state device opera¬ 
tion. The tracking of slow device dynamics is extremely 
challenging for the rigorous quantum kinetic approaches. 
Lastly, limiting the computational demand allows other 
details of an atomtronic circuit to be included, such as 
those involving atom injection and extraction processes, 
as well as complex confinement potentials. The end prod¬ 
uct is a predictive and flexible model that can be used for 
designing atomtronic devices and analyzing experimental 
data. 

Using the model, we address BEC phase-locking in a 
system that is simplified, but also geared more towards 
circuitry than the fully isolated system of Fig. 0 Figure [2] 
shows the model potential used for our current study. 
The system consists of two, weakly coupled harmonic 
wells of different frequency where atoms are injected into 
the left-hand well and extracted from the right-hand well. 
Additionally, atoms are extracted from the highest en¬ 
ergy states in each well analogous to the process of evapo¬ 
rative cooling used in typical BEC experiments. Notably, 
establishment of a phase relationship between coupled 
BECs has been studied previously in a similar geome¬ 
try m where dissipation and atomic interactions were 
found to play a key role in the phase-locking mechanism. 
In this work we study this phase evolution under the in¬ 


fluence of atom injection; and, as will be shown below, 
the strength of coupling between the ground states of 
each well strongly dictates the steady-state behaviour of 
the system. 

The remainder of the paper is organized as follows. 
Section II provides a description of the theory and model 
development. The working many-body equations of mo¬ 
tion are presented and the collision model described. Sec¬ 
tion III illustrates an application of the model to study 
atom transport between two adjacent harmonic poten¬ 
tials in the presence of atom injection and extraction. 
The simulations show two distinct dynamical scenarios 
involving the phase-locking of states in the different wells. 
Section IV summarizes the results and discusses future 
improvements to the model. 


II. THEORY 


The system Hamiltonian can be written as the sum of 
two parts: 


H = 


p ? 

t + VM 


+ )V{r i: rj), (1) 


i>j 


where pi is the momentum of the i th atom (all with mass 
to), V (r i) is the confining potential, and V (r i,rj) de¬ 
scribes interatomic interactions, which play an impor¬ 
tant role in the population dynamics of the system and 
are often described by a contact interaction mm- ai - 
though the model readily describes an arbitrary number 
of wells, for brevity we present the derivation assuming 
a double-well potential. In this case, we introduce the 
Hamiltonians for the uncoupled well potentials, V a ( r ) 
and Vb (r): 


Ha 


H b 


2m 
P 2 
2 TO 


+ V a (r), 
+ Vb (r), 


which have the solutions 


( 2 ) 


Ha4>? (r) = ef<j)f (r), 

H B 4>f (r) = ef 4>f (r) . (3) 

With V a (r) and 14 (r) given by the dashed parabolas de¬ 
picted in Fig. iEq- @ describes the familiar eigenvalue 
equation of the quantum harmonic oscillator. 

Using these eigenstates, we introduce the second- 
quantized wavefunction of the double-well system: 


$ ( r ) = ai ^ W + ( r ) > ( 4 ) 

i 3 

where (pf (r) and (pf (r) are localized in the left- and 
right-hand wells, respectively, and cq and bj are the an¬ 
nihilation operators for atoms in the i th and j th state 
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of the left- and right-hand wells, respectively. The cre¬ 
ation and annihilation operators obey the commutation 
relations, 


\cii ? ] — \pi ? bj ] 

= 0, 

(5) 

[ a i> bj} — 

a-i, &]■ 

= 0 , 

(6) 

a%,a ] = 

b h t 
°j 

= $i,j- 

(7) 


the following equations of motion for the system: 


lh it 

[4 a 4) 


- E 9ri - (ajbjj) , 

.7 

( 12 ) 

ih— 

dt 
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= E 9i:> (( b]at ) - (albjYj , 

i 

(13) 
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dt 
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Thus, in second-quantized form, 


A f 

a ik \ a i+k a i+k 


4bj 


H = [ d 3 r $ (r) -1 - V (r) ip (r) 

J L 2m J 

+ J d 3 r' J d 3 r ft (r') $ (r) V (r, r 7 ) ip (r) ip (r') , ( 8 ) 


giving essentially the two-site Bose-Hubbard Hamilto¬ 
nian [T 2 ] : 

H = E e 4 a l a i + E e ? b ) b i 

i j 

+ E {9ij a l b j + H.a4j 

ij 

+ E ( v m a i-k a ]+k a J a i + H - a ) 

i,j,k^ 0 

+ E (■ v fjk4-k b ]+k b j a i + H - a ■) 

+ E {<k b \-k b \ + k b jbi + H.a), (9) 

i,j,k^ 0 


~9ij (( b ] b j) - ( a !«*)) . ( 14 ) 

where Aj, = e® — ef is the energy shift between the left- 
and right-hand wells, and <j“ fc = v^ knk + v™ +k _ k , 
with a = A, B, is the coefficient for the energy renormal¬ 
ization. 

The final step in developing our model is to include 
the effects of collisions. Following a treatment that has 
been successful in describing collisional effects in semi¬ 
conductor devices, in terms of reproducing the results of 
quantum kinetic methods mm, we make an effective 
relaxation rate approximation and write 

= - 7 [< - , (! 5 ) 

col 

where 7 is an effective collision rate that drives the sys¬ 
tem towards an equilibrium distribution. Here, nf is 
the actual population of the i th level in the a well and 
is the equilibrium population that it relaxes 
to in the limit 7 —► 00 . In our case, /(ef,//,T) is the 
Bose-Einstein distribution function for a given chemical 
potential, /j,, and temperature, T, which are determined 
by number and energy conservation: 


<K_ 

dt 


where gtj describes a direct coupling between the two 
wells and Vij k accounts for interatomic correlations. 
These two coefficients are given by 


k +v v 


9ij = J d 3 r <pf (r) 

= J d 3 r' J d 3 r cpf_ k (r') <jjP. +k (r) 


( r ) i ( 10 ) 


a, 8 
V ijk 


xH (r, r') tpft (r) </>“ (r'), 


( 11 ) 


where a,fi = A, B. 


A. Equations of Motion and Collision Effects 


Using Eq. ([ 9 ]) and commutation relations for the 
atomic creation and annihilation operators, we arrive at 


E E<= E E/«^> r )> ( 16 ) 

a—A,B i <t—A,B i 

E E«= E ( 17 ) 


t=A,B 


7—A.B 


To picture the collision effects described by Eq. (151, 
consider again the scenario of Fig. [l|a). Assuming the 
system is completely isolated, the energy of the atoms 
can only be redistributed by interatomic collisions. Ini¬ 
tially, the atomic distribution is far from equilibrium as 
is evident from the right-hand plot in Fig. |lja) , where 
there are empty states in the population distribution of 
the two-well system at energies corresponding to levels in 
the right-hand well (e.g., at Em). As the system evolves, 
collisions drive the system towards equilibrium by redis¬ 
tributing the population among the available states. The 
end result is an equilibrium population distribution with 
higher lying states populated (i.e., at a higher effective 
temperature (see Fig. [ljb)) because of energy conserva¬ 
tion). This is the role of the collisions represented by 
Eqs. @ 
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III. APPLICATION OF THE MODEL 


To illustrate the application of the above model, we 
consider what might represent a portion of a device that 
is used in an atomtronic circuit. It consists of an asym¬ 
metric double-well [201 HU in which atom flux arrives at 
one well and leaves through the other. Of interest is the 
path by which this simple system arrives at a dynamical 
steady-state, while also remaining out of thermal equi¬ 
librium. 

For this example, a further simplification is made by 
using the following relations 


a] a, 


4 b i 


= n 


= 71a 


l 5 

B 


= yjnfnf ex p [-*%], 


(18) 

(19) 

( 20 ) 


where dij = Of — Of is the relative phase between states 
of the left- and right-hand wells. Inclusion of collision 
effects, as well as source and extraction contributions, P, 
and yj 4 ’" 8 , respectively, gives 


dnf 

dt 


dnl 

dt 


j CQl 

+Pi - ifnf, ( 21 ) 

i co1 

( 22 ) 


dOij 

dt 


9ij 


j t . 
fc^o 

(nf - nf) 


[pjk^j+k ®ik^i-\-k) 


nfnf 


cos 9 


■iji 


(23) 


In the simulations we consider the ground state of the 
left-hand well initially occupied together with a negligible 
background population in the other states (i.e., initially 
an essentially zero temperature BEC exists in the left- 
hand well). We consider six energy levels in the left- 
hand well and four in the right-hand well. The initial 
populations are nf = 0.08 and nf = 10~ 6 for i > 1 in 
the left-hand well, and nf = 0.001 in the right-hand well. 
We use an effective collision rate of 7 = 4 s -1 . Atoms 
are injected into the fifth level of the left-hand well at a 
rate of P 5 = 0.1 s -1 , and evaporated from the highest 
levels at a rate of = yf = 4 s -1 . The population in 
the lowest three energy levels of the right-hand well are 
outcoupled at a rate of y ® 2 3 = 0.1 s -1 . The temporal 
evolution of the populations in each well are studied in 
the presence of collisions and a direct coupling between 
the ground states of each well as described by Eq. (10) 
for 511 . 
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FIG. 3. Time evolution of the relative phase (solid, blue) 
and detuning (solid, red) between ground state populations 
of the two wells for the unlocked case. 


In this double-well system the presence of a relative 
phase between BECs in the left- and right-hand wells is 
of particular interest. Based on Eq. (23) phase-locking 
(i.e. d0n/dt = dtp/dt —> 0 ) can occur when the term 
containing gn in Eq. (23) is larger than the first term on 
the right-hand side of the same equation. The effect of 
phase-locking on this system is discussed in more detail 
below where the time evolution of the system is studied 
for different coupling strengths between the two wells of 
the potential. 


A. Unlocked operation 


We begin by considering a ground state coupling pa¬ 
rameter of gn/h = 1 s -1 . Figure [ 3 ] shows the temporal 
evolution of ip. Also plotted is the ground state frequency 
difference, or detuning, which is given by 



A n + (vfk n i+k - vik n t+k) . 
fc^O 


(24) 


where crf k /h = crf k /h = 0.2 s 1 for all k 0. For 
911/ft = 1 s-\ the BEC phases are unlocked and the 
relative phase between the two wells exhibits rapid oscil¬ 
lations. Additionally, the detuning remains nearly con¬ 
stant in time. 

Figure [4] shows the time evolution of the spatial dis¬ 
tribution of the ensemble-averaged atomic population, 
which is given by 



+2\J~nfnf Re {<j>f<j>f exp [ip]] . (25) 

The initial population distribution is localized in the 
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FIG. 4. Spatial distribution of the ensemble-averaged atomic 
population, (ip'ip}, in the unlocked case for t = 10~ 6 s (solid, 
black), 1 s (dotted, blue), 624.5 s (dash-dotted, green), and 
626.4 s (dashed, red). At large times the population distri¬ 
bution undergoes small oscillations between the two wells as 
depicted by the dash-dotted (green) and dashed (red) curves. 


left-hand well as depicted by the t = 10~ 6 s curve. The 
remaining curves show the temporal evolution of the pop¬ 
ulation distribution where small oscillations between the 
non-degenerate eigenstates <p± and (pf occur, as depicted 
by the t = 624.5 s and 626.4 s curves. 


B. Locked operation 

System behavior becomes distinctly different when 
gu/h increases to 2 s -1 . Figure [S] explicitly shows the 
phase-locking, which occurs after a time t > 2 s. The 
locking arises from an increase in the ground state pop¬ 
ulation, n) 4 , of the left-hand well due to the pump and 
extraction processes. This in turn causes the exchange 
shift to decrease the detuning between left- and right- 
hand well ground states (red curve, Fig. [5|. When the 
detuning becomes sufficiently small, phase locking occurs 
(compare blue curves, Figs. [3] and [5]). 

Figure [6] shows the evolution of the same initial atomic 
spatial distribution as in the unlocked case. Here, the 
populations in both wells continue to grow during the 
transfer of atoms from the left-hand to the right-hand 
well, until a time-independent distribution is reached, 
with substantially higher populations in both wells than 
in the unlocked situation (compare to Fig. [4]). Underly¬ 
ing the time-independence and population growth is the 
formation of a composite eigenstate, resulting from the 
phase locking of the ground states in both wells. 



Time (s) 

FIG. 5. Time evolution of the relative phase (solid, blue) 
and detuning (solid, red) between ground state populations 
of the two wells for the locked case. 



z (pm) 

FIG. 6. Spatial distribution of the ensemble-averaged atomic 
population, (ip'ip \, in the locked case for t = 10“ 6 s (solid, 
black), 1 s (dotted, blue), 50 s (dash-dotted, green), and 
1000 s (dashed, red). A time-independent distribution is 
reached where an appreciable population has built up in the 
left well. 


C. Comparison of locked and unlocked operation 

The significance between locked and unlocked opera¬ 
tion for circuit operation is clearly evident when observ¬ 
ing the number of atoms leaving the system. In Fig. [7] 
the atom flux leaving the right-hand well is plotted as a 
function of time. When the system is not phase-locked 
the output flux fluctuates around an average value of 
~ 0.04 atoms/s. The fluctuations arise as a result of 
a small back and forth population transfer between the 
two wells. When phase-locking occurs, the output flux 
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FIG. 7. Time evolution of the output flux for the unlocked 
(solid, blue) and locked (solid, red) case. 



FIG. 8. Steady-state population distributions in the un¬ 
locked (solid, blue) and locked (solid, red) cases. Each data 
point corresponds to the population at a given energy level. 
The solid lines show the Bose-Einstein distributions with 
equivalent total population and energy. 


increases by a factor of 3.5 to a steady-state value of 
~ 0.14 atoms/s. 

Some insight into the increase in output flux due to 
phase-locking can be gained by studying the population 
in each level once steady-state is reached. Figure [8] shows 
the population distribution in steady-state for both the 
unlocked and locked cases. As time progresses in the 
unlocked case, relaxation takes place and the popula¬ 
tions settle essentially to a Bose-Einstein distribution af¬ 
ter t ~ 10 s. Note that the population distribution is 
not exactly Bose-Einstein in nature due to the effects 
of injecting and extracting atoms at specific energy lev¬ 
els. A similar progression is followed at early times in 




FIG. 9. Time evolution of the (a) fractional ground state 
population and (b) effective temperature of atoms in the left- 
hand well for the unlocked (solid, blue) and locked (solid, red) 
cases. 


the locked case; however, with the onset of phase-locking 
substantial cooling takes place (indicated by the decrease 
in population of higher lying states) while the total atom 
population continues to grow appreciably. Time inde¬ 
pendence in this case occurs around t ~ 100 s. 

Fig. 0 summarizes the results presented in Fig. [H by 
showing the time evolution of the fractional ground state 
population and effective temperature of atoms in the left- 
hand well. Here, the fractional ground state population 
in the left-hand well is given by nj 4 / n t • This quan¬ 

tity is plotted in Fig. |9](a), where the difference in con¬ 
densate occupation is clearly depicted. Figure[9](b) shows 
the effective temperature of the left-hand well, which is 
obtained by fitting the population distribution to a Bose- 
Einstein distribution. At early times, both locked and 
















unlocked systems show a rise in effective temperature as 
chemical potential is converted to thermal energy and 
the system trends towards equilibrium. In the unlocked 
case the system appears to reach a steady-state with no¬ 
ticeable residual fluctuations. On the other hand, in the 
locked case the system moves towards a time-independent 
operation exhibiting significant cooling and BEC growth 
in the left-hand well. This growth leads to the large 
increase in the number of atoms flowing out of the right- 
hand well as shown in Fig. [7] 

The physics of coherent and incoherent tunneling has 
drawn interest since the 1920’s. In regards to macro¬ 
scopic quantum tunneling experiments with BEC, recent 
reports have also indicated two regimes of operation: 
the Josephson oscillation regime and the self-trapping 
regime [221 . which are determined by the relative phase 
of two BECs occupying neighboring potential wells. Di¬ 
rect relation of our modeling results to experimental work 
is ongoing. The connection is challenging because our 
model casts system behavior in terms of the energy eigen¬ 
states of the wells and how they are effected by popula¬ 
tion dynamics, such as the exchange shifts. In contrast, 
the experimental work relates system behavior directly 
to the atomic populations. 

IV. CONCLUSION 

A model is developed for use in the design of 
atomtronic circuitry based on finite temperature Bose- 
condensed gases. Working in the Heisenberg picture, 
equations of motion for atomic populations and coher¬ 
ences are derived. Energy renormalizations are treated 
at the mean-field level, and collision effects are taken 
into account using an effective relaxation rate descrip¬ 
tion. Numerical solution of the population and coherence 
equations of motion allows tracking of the dynamics of fi¬ 
nite temperature BECs in the presence of atom injection 


and extraction. 

The model is demonstrated by studying the evolution 
of a Bose-condensed gas in the presence of atom injec¬ 
tion and extraction in a double-well potential. In the 
presence of collisions and dissipation, two regimes of de¬ 
vice operation are observed corresponding to unlocked 
and locked operation. These operating regimes are iden¬ 
tified based on whether the population transfer between 
the two wells of the potential is incoherent (unlocked) 
or coherent (locked). Behavior of the system in terms 
of the output flux, degree of condensation, and effective 
temperature is shown to be distinctly different for each 
operating regime. 

The goal of this paper is to introduce a predictive 
and flexible model for use in designing atomtronic de¬ 
vices and for analyzing BEC experiments. There is much 
room for further improvement of the model. For exam¬ 
ple, quantum fluctuations that are especially important 
at very low temperature may be included by extending 
the equations of motion to include higher order corre¬ 
lations described by terms like (a\a\aia^j, (bjbjbjbj\, 

and (^a\bjbjai^. In addition, the role of dissipation and 
dephasing may be treated more consistently by use of the 
equivalence to the Lindblad terms in the density operator 
approach [[25] , 
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